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In the framework of QCD sum rules, we present an improved study of our previous work [Phys. 
Rev. D 80, 056004 (2009)] particularly on the DD* molecular state to investigate that the possibility 
of the newly observed Z c (3900) as a 5- wave DD* molecular state. To ensure the quality of QCD sum 
rule analysis, contributions of up to dimension nine are calculated to test the convergence of operator 
product expansion (OPE). We find that the two-quark condensate (qq) is very large and makes the 
standard OPE convergence (i.e. the perturbative at least larger than each condensate contribution) 
happen at very large values of Borel parameters. By releasing the rigid OPE convergence criterion, 
one could find that the OPE convergence is still under control. We arrive at the numerical result 
3.86±0.27 GeV for DD* , which agrees with the mass of Z c (3900) and could support the explanation 
of Z c (3900) in terms of a S-wave DD* molecular state. 
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I. INTRODUCTION 

Very recently, BESIII Collaboration studied the process e + e — > tt + tt J/ip at a center-of-mass energy 
of 4.26 GeV, and reported the observation of a new charged charmonium-like structure Z c (3900) in the 
^ x ± J/ip invariant spectrum with a mass of 3899.0±3.6±4.9 MeV and a width of 46±10±20 MeV [l|. Before 
the BESIII's observation, Chen et ai. predicted that a charged charmonium-like structure is observable in 
the V(4260) ->• J/ipir+ir- process [2]. In the study of V(4260) -> tt+tt" J/ip decays, Belle Collaboration 
also observed a £(3895)+ state with a mass of 3894.5 ± 6.6 ± 4.5 MeV and a width of 63 ± 24 ± 26 MeV 
in the 7r ± J/ip mass spectrum [3j. Xiao et al. confirmed the charged state Z c (3900) in the analysis of data 
taken with the CLEO-c detector at -0(4160), and measured its mass and width to be 3885 ± 5 ± 1 MeV 
and 34 ± 12 ± 4 MeV, respectively 0|. 

The new experimental results have aroused theorists' great interest in comprehending the internal struc- 
tures of Z c (3900). Soon after the observation of Z c (3900), it was proposed that these states are S-w&ve 
DD* molecules [5J, l6j ■ Subsequently, there also appeared many other works to explain this exotic states 
j7Hl3j. Undoubtedly, it is interesting and significative to investigate that whether Z c (3900) could be a 
S- wave DD* state. To understand the inner structure of i? c (3900), it is very helpful and quite needed to 
determine their properties like masses quantitatively. Nowadays, QCD is widely believed to be the true 
theory of describing strong interactions. However, it is quite difficult to acquire the hadron spectrum from 
QCD first principles. The main reason is that low energy QCD involves a regime where it is futile to 
attempt perturbative calculations and the strong interaction dynamics of hadronic systems is governed by 
nonperturbative QCD effects completely. Meanwhile, one has limited knowledge on nonperturbative QCD 
aspects for there are still many questions remain unanswered or realized only at a qualitative level. 



The method of QCD sum rules [14j is a nonperturbative formulation firmly based on the basic theory of 
QCD, which has been successfully applied to conventional hadrons (for reviews see 15l4l8| and references 



therein) and multiquark states (e.g. see [19j). In particular for the S-w&ve DD* molecular state, we have 



definitely predicted its mass to be 3.88 ± 0.10 GeV with QCD sum rules several years ago in Ref. [20j, in 
which mass spectra of molecular states with various {Qq\{Q^ 'q} configurations have been systematically 
studied. Numerically, one could see that our prediction for the mass of DD* state agrees well with the 
experimental data of newly observed Z c (3900). That result could support the explanation of Z c (3900) as a 
S- wave DD* molecular state. At present, we would put forward an improved study of our previous work on 



the DD* state in view of below reasons. First, it is known that one can analyze the OPE convergence and 
the pole contribution dominance to determine the conventional Borel window in the standard QCD sum 
rule approach to ensure the validity of QCD sum rule analysis. However, we find that it may be difficult to 
find a conventional work window rigidly satisfying both of the two rules in some recent works 21], which 



actually has also been discussed in some other works (e.g., Refs. (22H24I ]]). The main reason is that some 
high dimension condensates are very large and play an important role in the OPE side, which makes the 
standard OPE convergence happen only at very large values of Borel parameters. Whereas in the previous 
Ref. [20J , we merely considered contributions of the operators up to dimension six in OPE and the Borel 
windows are roughly taken the same values for the similar class of states for simplicity and convenience. 
Thus, it may be more reliable to test the OPE convergence by including higher dimension condensate 
contributions than six and considering the work windows minutely, and then one could more safely extract 
the hadronic information from QCD sum rules. Second, even higher condensate contributions may not 
radically influence the character of OPE convergence in some case, one still could attempt to improve the 
theoretical result because some higher condensates are helpful to stabilize the Borel curves. Particularly 
for the newly observed Z c (3900) states, they can not be simple cc conventional mesons since they are 
electric charged. It may be a new hint for the existence of exotic hadrons and Z c (3900) are some ideal 
candidates for them. Once exotic states can be confirmed by experiment, QCD will be further tested and 
then one will understand QCD low-energy behaviors more deeply. Therefore, it is of importance and worth 
to make meticulous theoretical efforts to reveal the underlying structures of Z c (3900). All in all, we would 
like to improve our previous work to investigate that whether Z c (3900) could serve as a DD* molecular 
state. 

The rest of the paper is organized as follows. In Sec. [TT1 QCD sum rules for the molecular states are 
introduced, and both the phenomenological representation and QCD side are derived, followed by the 



numerical analysis and some discussions in Sec. Mil The last Section is a brief summary. 

II. MOLECULAR STATE QCD SUM RULES 

The starting point of the QCD sum rule method is to construct a proper interpolating current to represent 
the studied state. One knows that the method of QCD sum rules has been widely applied to multiquark 
systems since the experimental observations of many new hadrons in these years. At present, currents 
of molecular states and tetraquark states could be differentiated by their different construction ways. 
Concretely, molecular currents are built up with the color-singlet currents of their composed hadrons 
to form hadron-hadron configurations of fields, which arc different from currents of tetraquark states 
constructed by diquark-antidiquark configurations of fields. What's needed to note is that these two types 
of currents can be related to each other by Fiertz rearrangement. However, the transformation relations 
are suppressed by corresponding color and Dirac factors [19( and one could obtain a reliable sum rule while 
choosing the appropriate current to represent the physical state. This means that if the physical state is 
a molecular state, it would be best to choose a meson-meson type of current so that it has a large overlap 
with the physical state. Similarly for a tetraquark state, it would be best to choose a diquark-antidiquark 
type of current. When the sum rule reproduces a mass consistent well with the physical value, one can 
infer that the physical state has a structure well represented by the chosen current. In this way, one can 
indirectly and commonly discriminate between the molecular and the tetraquark structures of observed 
states. One can expect that the judgements could be very effective for some ideal cases, e.g. the results 
obtained from different types of currents are very different so that one could easily discriminate them. Note 
that in some exceptional cases that it may not be very different for final results from molecular currents 
and tetraquark currents. For example, Narison et al. investigated both molecular and tetraquark currents 
associated with X(3872) and they finally gained the same mass predictions within the accuracy of QCD 
sum rule method in Ref. [25]. For the present work, in order to study the possibility of Z c (3900) as a 



S'-wave DD* molecular state, we thus construct the molecular current from corresponding currents of D 
and D* mesons to form meson-meson configurations of fields. In full theory, the interpolating currents for 
heavy D mesons can be found in Ref. [261 ]. Therefore, one can build the following form of current 

3% D * = {Qaii5qa)(qb'i ll Qb), 

for DD* with J p = 1 + , where q indicates the light u or d quark, Q denotes the heavy c quark, and the 
subscript a and b are color indices. Note that the quantum numbers of Z c (3900) have not been given 
experimentally for the moment, and 1 + is just one possible choice of their spin-parities. 
One can then write down the two-point correlator 

IT"V) =iJ<Px***(0\T\j& D .(x)j u B + D .(0)]\0)- (!) 

Lorentz covariance implies that the correlator can be generally parameterized as 

n^V) = (^-^V^ 2 ) + ^n(°V). (2) 

The term proportional to g^ will be chosen to extract the mass sum rule. Phenomenologically, IT 1 '^ 2 ) 
can be expressed as 

TTftV 2^ [A (1) ] 2 1 r , ImnWP 11011 ^) 

n (1) (<T) = —^ — n + - / ds ^— — + subtractions, (3 

M DD' 'I * J so $ - 1 

where Mjj D * denotes the mass of the DD* state, sq is the continuum threshold parameter, and \^ gives 
the coupling of the current to the hadron (0|j^„, \DD*) = A^e^. In the OPE side, H^(q 2 ) can be written 

as 

nU(q 2 )= ds p _ \ s) +rrf m V), (4) 

J4m^ s 1 

where the spectral density is given by p (s) = -ImllA-v(s). Technically, we work at leading order in 



a s and consider condensates up to dimension nine, employing the similar techniques as Refs. 27|, l28l |. 
To keep the heavy-quark mass finite, one can use the momentum-space expression for the heavy-quark 
propagator [26] 

S Q(P) = 1 t.9^G^ a (0)-^ j-2 [ct kA (^ + m Q ) + {i> + m Q )a K \] 



- l -g 2 t A t B Gip(0)GfM I ± m l [ 7 °(^ + m Q h^ + mqWU + m Q )Y (5) 

4 {J> m Q) 

+ l a {i> + m Q)l^(^> + ™q)7 /3 (^ + WqIt" + l a ii> + to q)7 M (^ + m Q)l U {i> + m Qh^](p' + m Q) 
+ ^g z f ABC G^GlG% {p2 _ X m2 $ + m Q Mp 2 - 3m 2 Q ) + 2m Q (2p 2 - m%)](j + ™q)- 

The light-quark part of the correlator can be calculated in the coordinate space, with the light-quark 
propagator 

T O h ?T iJ J. TO h 

- Yi^igqv ■ Gq) + 27 _° 2 m q (gqa ■ Gq)fi ~ ^ ^ 3 (qq)(g 2 G 2 ), (6) 

which is then Fourier-transformed to the momentum space in D dimension. Since masses of light u and d 
quarks are three order magnitudes less than the one of heavy c quark, they are neglected here following 



the usual treatment. The resulting light-quark part is combined with the heavy-quark part before it is 
dimensionally regularized at D = 4. After equating Eqs. ([3]) and Q, assuming quark-hadron duality, and 
making a Borel transform, the sum rule can be written as 



[A (i) ]2e -M? D ,/A/ 2 = f S0 dsp OVV e -s/M> + ^ncond^ 



(7) 



with M 2 the Borel parameter. Making the derivative in terms of M 2 to the sum rule and then dividing 
by itself, we have the mass of DD* state 
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where 



p 0PE (s) = pP ert ( S ) + p<™>(«) + ptf°')(8) + p^- G ^(s) + p^\s) + p^ G3 \s) + p<^ 2G2 >( s ), 



with pP ert , p^\ p(9 2c2 \ p (9<i<r-Gq) ) p (qq)\ p (a 3 G 3 )^ and p {qq){g 2 G 2 ) are the per turbative, two-quark con- 
densate, two-gluon condensate, mixed condensate, four-quark condensate, three-gluon condensate, and 
two-quark multiply two-gluon condensate spectral densities, respectively. In fact, the spectral densities up 
to dimension six have been given in our previous work [20], which are also enclosed here for the paper's 
completeness. Concretely, the spectral densities are 
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with (qq)(gqa ■ Gq), (g 2 G 2 ) 2 , (qq)(g 3 G 3 ), and (g 2 G 2 )(gqa ■ Gq) denoting the two-quark multiply mixed 
condensate, four-gluon condensate, two-quark multiply three-gluon condensate, and two-gluon multiply 
mixed condensate, respectively. 

III. NUMERICAL ANALYSIS AND DISCUSSIONS 

In this section, the sum rule ([5]) will be numerically analyzed. The input values are taken as m c = 1.23 ± 
0.05 GeV, (qq) = -(0.23 ± 0.03) 3 GeV 3 , (gqa ■ Gq) = m 2 (qq), mg=0.8±0.1 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 , 



and (g 3 G 3 ) — 0.045 GeV 6 16|. In the standard procedure of sum rule analysis, one should analyze the OPE 
convergence and the pole contribution dominance to determine the conventional Borel window for M 2 : on 
the one side, the lower constraint for M 2 is obtained by considering that the perturbative contribution 
should be larger than each condensate contribution to have a good convergence in the OPE side; on the 
other side, the upper bound for M 2 is obtained by the consideration that the pole contribution should be 
larger than the continuum state contributions. At the same time, the threshold ^/sq is not arbitrary but 
characterizes the beginning of continuum states. Hence, the most expected case is that one could naturally 
find the conventional Borel windows for studied states to make QCD sum rules work well. 

In order to test the convergence of OPE, its various contributions, i.e. the perturbative, two-quark, 
two-gluon, mixed, four-quark, three-gluon, two-quark multiply two-gluon, two-quark multiply mixed, four- 
gluon, two-quark multiply three-gluon, and two-gluon multiply mixed condensate contributions, are com- 
pared as a function of M 2 and showed in FIG. 1. Graphically, one could see that in the OPE side there 



exists some similar problem which has been discussed in some of our recent works 2l| and others e.g. 



22 



24j |. Concretely, here some condensates especially two-quark condensate (qq) are very large and play 
an important role in the OPE side, which makes the standard OPE convergence (i.e. the perturbative at 
least larger than each condensate contribution) happen only at very large values of M 2 . The consequence 
is that it is difficult to find a conventional Borel window where both the pole dominates over the continuum 



and the OPE converges well. Following the similar treatment in Refs. |21|, we could try releasing the rigid 
convergence criterion of the perturbative contribution larger than each condensate contribution here. It is 
not too bad for the present case, there are two main condensates i.e. (qq) and (gqa ■ Gq) and they could 
cancel out each other to some extent since they have different signs. What is also very important, most of 
other condensates calculated are very small, which means that they could not radically influence the char- 
acter of OPE convergence. Therefore, one could find that the OPE convergence is still under control here. 
In the phenomenological side, the comparison between pole and continuum contributions of sum rule ([7]) 
as a function of the Borel parameter M 2 for the threshold value ^/Iq — 4.4 GeV is shown in FIG. 2, which 
shows that the relative pole contribution is approximate to 50% at M 2 = 2.7 GeV 2 and decreases with 
M 2 . Similarly, the upper bound values of Borel parameters are M 2 = 2.6 GeV for y/so — 4.3 GeV and 
M 2 = 2.9 GeV 2 for </sE = 4.5 GeV. Thus, the Borel window for DD* is taken as M 2 = 2.1 - 2.7 GeV 2 
for y/so~ = 4.4 GeV. Similarly, the proper range of M 2 are 2.1 ~ 2.6 GeV 2 for y/s^ = 4.3 GeV and 
2.1 ~ 2.9 GeV 2 for ^/sq = 4.5 GeV. The mass of the DD* molecular state as a function of M 2 from sum 
rule {8} is shown in FIG. 3 and it is numerically calculated to be 3.86 ±0.13 GeV in the above chosen work 
windows. Considering the uncertainty rooting in the variation of quark masses and condensates, we gain 
3.86 ± 0.13 ± 0.14 GeV (the first error reflects the uncertainty due to variation of ^/s ~ and M 2 , and the 
second error resulted from the variation of QCD parameters) or concisely 3.86 ± 0.27 GeV for the S-wave 
DD*. 
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FIG. 1: The OPE contribution in sum rule ([Tj) for ^/io = 4.4 GeV. The OPE convergence is shown by comparing 
the perturbative, two-quark condensate, two-gluon condensate, mixed condensate, four-quark condensate, three- 
gluon condensate, two-quark multiply two-gluon condensate, two-quark multiply mixed condensate, four-gluon 
condensate, two-quark multiply three-gluon condensate, mixed multiply two-gluon condensate contributions. 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 

°r 









— pole/total 

— continuum/total 


o 

CO 

M 2 (GeV 2 ) 


,''' 



2.5 



3.5 



FIG. 2: The phenomenological contribution in sum rule Q for t/so — 4.4 GeV. The solid line is the relative pole 
contribution (the pole contribution divided by the total, pole plus continuum contribution) as a function of M 2 
and the dashed line is the relative continuum contribution. 



IV. SUMMARY 



Stimulated by the newly observed charged charmonium-like structure Z c (3900) for which can not be 
simple cc conventional mesons and are some ideal candidates for exotic hadrons, we present an improved 
QCD sum rule study of our previous work on the DD* molecular state to investigate that whether it could 
be a S-wave DD* molecular state. In order to ensure the quality of QCD sum rule analysis, contributions 
of up to dimension nine are calculated to test the convergence of OPE. We find that some condensates 
in particular (qq) play an important role and make the standard OPE convergence (i.e. the perturbative 
at least larger than each condensate contribution) happen at very large values of Borel parameters M 2 . 
By releasing the rigid OPE convergence criterion, one could find that the OPE convergence is still under 
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FIG. 3: The mass of the DD* molecular state as a function of M 2 from sum rule (0. The continuum thresholds 
are taken as y/sE = 4.3 ~ 4.5 GeV. The ranges of M 2 is 2.1 ~ 2.6 GeV 2 for y/s^ = 4.3 GeV, 2.1 ~ 2.7 GeV 2 for 
^ = A A GeV, and 2.1 ~ 2.9 GeV 2 for y/so = 4.5 GeV. 



control and the final result 3.86 ± 0.27 GeV is obtained for the S-wave DD* molecular state, which 
coincides with the experimental data of Z c (3900). From the final result, one could assuredly state that it 
could provide some support to the DD* molecular explanation of Z c (3900). At the same time, one should 
note that the DD* molecular state is just one possible theoretical interpretation of Z c (3900), and it does 
not mean that one could arbitrarily excluded some other possible explanations (e.g. tetraquark states) at 
the present time just from the result here. In fact, more minute information on the nature structures of 
Z c (3900) could be revealed by the future contributions of both experimental observations and theoretical 
studies. 

Note added- As we were preparing to submit this paper, we became aware of a paper from our colleague 
that also analyzes Z c (3900) as a DD* molecular state with QCD sum rules [29j, but then they consider 
contributions up to the same dimension six as our previous work [20(. 
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